Destruction of Peierls dimerization in the molecular crystal model: 
Effects of quantum phonon fluctuations 
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Effects of quantum pfionon fluctuations on the Peierls dimerization in the one-dimensional molecular 
crystal model are reexamined by a functional integral approach. An equation for the dimerization 
order parameter is obtained within a one- loop approximation. The destruction of Peierls dimer- 
ization is found by taking the quantum phonon fluctuations into account for the spin- 1/2 electron 
system. The consistency and discrepancy of the result with previous ones are discussed. 
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In a one-dimensional electron-phonon system, the 
metallic phase is unstable against the electron-phonon 
interaction, which results in the Peierls transitionE At 
the half-filling dimerized lattce will be formed 

and the ground state is of a dimerized long-range or- 
der. This is true for an arbitrarily small electron- 
phonon coupling within the mean-field adiabatic approx- 
imation which treats the phonon degree of freedom clas- 
sically. It is an interesting problem whether the dimer- 
ized ground state survives the quantum phonon fluctua- 
tions. There have been many works to study this prob- 
lem. Among them are the renormalization groun-.anal- 
ysis plus the Monte Carlo numerical |^alculationaj'0 on 
the Su-Schrieffer-IIeegeij_(SSH) modelel and the molecu- 
lar crystal (MC) model,El the variational approachQu on 
the same niDdels, and the one-loop approximation on the 
SSH modelEI and the halogen- bridged mixed-valence lin- 
ear chain complexes (MX) modeLLi 

In the SSH model as well as the MX model, the op- 
tical phonon-electron interactions are dominant and the 
acoustic phonon-electron interactions are weak, so the 
quantum fluctuations are mainly induced by the opti- 
cal phonon-electron couplings. It has been shown that 
the dimerization of lattice survives the quantum phonon 
fluctuations for the spin^i systems whereas it is broken 
for the spinless systems.HiJ However, in the MC model, 
the interactions between electrons and two branches of 
phonons are equally important, it seems that the quan- 
tum fluctuations should be much stronger than that in 
the SSH and MX models. Although the renormalization- 
group analysis shows absence of dimerized long range or- 
der for the zero mass system, Hirsch and Fradkin con- 
cluded that the ground state is dimerized for any nonzero 
mass systems from the analysis of symmetries as well as 
their Monte Carlo numerical simulations .3 

In this paper, we reexamine the effects of quan- 
tum phonon fluctuations on Peierls dimerizaton in the 
molecular-crystal model by a functional integral ap- 
proach. An equation for the dimerization order param- 
eter is obtained within a one-loop approximation. We 
find the Peierls dimerization will be destroyed by tak- 



ing the quantum phonon fluctuations into account for 
the spin-i electron system. An equivalency between the 
spin-i MC model and the spinless SSH or MX model is 
shown. The consistency and discrepancy of our results 
with the renormalization- group analysis as well as the 
Monte Carlo numerical simulations are discussed. 

The one-dimensional molecular-crystal madel we con- 
sider in this paper takes the following form.B 
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where c\ ^ and q^s are the creation and annihilation op- 
erators of electrons sitting at molecular / with spin s, ui 
(its conjugated momentum is pi) describes an internal 
degree of freedom of molecule I, to, a and K are the con- 
stants for the electron hopping, electron-phonon coupling 
and lattice elasticity, and M is the mass of the internal 
vibrations in molecules which form the one-dimensional 
lattice. The model can also be viewed as describing the 
coupling between electrons and excitons with boson char- 
acters. 

To take into account the quantum phonon fluctua- 
tions, we employ a functional integral approach. The 
partition function of the system can be written as a func- 
tional integral over both complex and Grassmann vari- 
ables: 
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where the action S is 
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By making Fourier transformations on the integral vari- 
ables 
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and similarly for (jj^ ^ and (f>i^s, then for the phonon vari- 
ables we have two branches on a dimerized lattice, that 
is, the acoustic phonon ag(f2)[= Uq{fl)] and the optical 
phonon bq{il)[= Uq+Q{fl)], Q = 2kp {kp = 7r/2 for the 
half-filled band case, the lattice constant has been set as 
unit). For the Grassmann variables we could introduce a 
spinor notation viaE3 
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which are right- and left-moving electrons, respectively, 
then we have the partition function as 
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and the corresponding action could be written as a sum 
of five terms 



S — Sa + iSfc + S'e + Sa- 
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where the acoustic- and optical-phonon parts of the ac- 
tion are, respectively. 
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Sb^Y. + W^''] bgmb_,{^n), (8b) 
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the electronic part of the action is 

Se^Y. *M(^)[(«w-A') + 2tosinfccj3]$fe,,(w), (8c) 

k,u,s 

and the acoustic and optical phonon-electron interaction 
parts are the following 
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Sb^, = E E bq{mijLu)ai<i>k-qA^ - (8e) 

where ai are Pauli matrices. 

In the SSH model and the MX model, the opti- 
cal phonon-electron interactions are dominant and pthe. 
acoustic phonon-electron interactions can be ignoredBt2l 
However, as seen from the Eq. (8d) and (8e), the cou- 
plings between electrons and two branches of phonons 
are equivalent in the MC model. The fact implies we 



should take the both into our consideration. One can in- 
tegrate out the electronic variables of the partition since 
the action is bilinear in the Grassmann fields, the result- 
ing partition function is 



Z = V[a*ab*b]e' 
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and the effective action is 



Scff = E [5-^ + hM^^] [aq^a^qi-^) + b,{n)b^,i-n)] 

q,n 

-iVlndet(Al), (10) 

where N is the spin degree of freedom, det denotes the 
determinant of the matrix M, which is defined by the 
actions Se + Sa-e + Sb-e — J2s ^IM^s and has fc, uj, 
and the indices of the Pauli matrices as labels. 

By the functional derivation of the effective action 
Scff with the phonon variables vanishing, we have the 
Peierls dimerization, i.e., the optical phonon condensa- 
tion at zero momentum and zero frequency {abq{Q)) = 
\/PLAdqflSn.Oi where the phonon order parameter A de- 
scribes the Peierls dimerization. Defining the noninter- 
acting electronic Green's function Go{k,Lu) by 

Go(fc, u}) — —{iuo — /i + 2io sinfc(T3 + Acri)^"^, (11) 

and bqiVl) = {bq{il)) + bq{fl), we obtain the effective ac- 
tion by expanding the logarithm in order of b and a: 
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is the zeroth-order contribution in b and a and is 



given by 
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in the thermodynamic limit, the dimensionless electron- 
phonon coupling constant A is defined by A = 
Na^ /nKto, the chemical potential ^ = at the 
half -filling case. The electronic spectrum Ek — 
±a/ (2io sin fc)^ + A^, the electronic gap is 2 A. 
The mean field gap equation 
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follows immediately from dS^^ /dA — 0, where A is the 
integral cutoff. The mean-field result tell us that the 
ground state is dimerized for any electron-phonon inter- 
actions. To investigate the effects of quantum phonon 
fluctuations, we should include higher order contribu- 
tions. The first-order term S^^ vanishes since S^g is ob- 
tained from the saddle point approximation. The second- 
order contribution to the effective phonon action is 
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where the functions f{q,fl) and g{q,fl) are defined as 
fohows 
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By performing the integration in Eq. (9) over the fluc- 
tuations b and a to the second-order term, we have 
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where the free energy density r[A] is composed of two 
parts, 



r[A] = ro[A] + ri[A], 



(18) 



To [A] — sf^ / PL is the zeroth-order contribution, and 
the one-loop contribution Fi [A] is given by 

ri[A] = i| ^^{x^[iMn' + xKf{q,n)] 

+ In [\A'm^ + XKg{q, f7)] } , (19) 

in the zero-temperature and thermodynamic hmit. The 
equation for the phonon order parameter A is determined 
by 9r[A]/i9A = 0. Within the one-loop approximation, 
we have 
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with both u = f and g, the renormalized phonon fre- 
quency fig = 2Au;q {ljq = y/K/M). In the limit of 
weak electron-phonon interaction, only electrons near the 
Fermi surfaces be important, that is, we can take small 
momentum approximation in the calculation of the func- 
tions /(g, ri) and g{q,V,), which gives f{q,i^) ~ g{q,il). 
This result implies the equivalence of the contributions 
from the optical and acoustic phonons in quantum fluc- 
tuations in the MC model. The derivation of the func- 
tion u{q, n) with the phonon order parameter A can be 



calculated straightforwardly, it gives u'{q,rt)/u{q,il,) 



Ru{q, n) = {(Pn^ + Atlq^ + 4A2)/4A, 



(22) 



with the abbreviation = 1 -|-4A2/riQ, the higher-order 
terms have been neglected as usual (in a renormalization- 
group study). The equation (20) can be compared with 
that in the SSH or MX model where only one Ru ap- 
pears since the acoustic phonon-electron interaction can 
be ignored, so that the spin-i {N — 2) MC model is 
equivalent to the spinless [N — 1) SSH or MX model in 
this sense. 

Now we obtain the phonon order parameter 
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by performing the integral of Eq. (20). It can be seen 
that in the adiabatic limit Af ^ oo, the quantum fluc- 
tuation is completely suppressed and the above equation 
returns to the mean-field one (14). 

The results obtained from Eq. (23) are discussed as 
follows. (1) For the spinless electron systems, the ground 
state is undimerized for a weak electron-phonon cou- 
pling whereas it is dimerized for a strong electron-phonon 
coupling. The numerical results are in goad agreement 
with that of the Monte Carlo simulationsE (2) For the 
spin-i electron systems, the ground state is undimerized 
in the limit of ill = 0; it is dimerized for the system 
with any nonzero mass in the limit of strong electron- 
phonon interactions. This result agrees with that of the 
renormalization-group analysis.a (3) A typical curve of 
our numerical results for the N ~ 2 case is given in Fig. 1. 
To compare the data we obtain with that of the Monte 
Carlo calculation, we used the phonon-staggered order 
parameter nip, which is related to A through nip = A/a. 
It is found that the Pcierls dimerization is destroyed in 
weak electron-phonon coupling regime. This is consis- 
tent with the spinless SSH or MX model since they are 
equivalent. The data we obtained are quite close to that 
of Monte Carlo simulation,El in which a jump in nip is 
seen at a critical a. We believe that this is a signa- 
ture of the transition discussed in the present work, al- 
though an opsite conclusion was reached at by the origi- 
nal authors. This result also implies that the break of the 
continuous symmetry in spin space does not necessarily 
produce dimerization although the continuous symmetry 
does prevent the existence of long- range dimerization 
order in the model.0 E'ipally, we mention the works by a 
variational approachjju which can not give out a transi- 
tion since it underestimates the quamtun fluctuations. 

In summary, the effects of quantum phonon fluctu- 
ations on Peierls dimerization in the molecular crystal 
model are reexamined by a functional integral approach. 
The calculation is performed by first integrating out the 
electronic variables of the partition function and then 
expanding the effective action to the quadratic terms in 
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the phonon variables. It is found that the Peierls dimer- 
ization is destroyed by the quantum phonon fluctuations 



superconductivity in this kind of systems. 
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FIG. 1. The dependence of the phonon-staggered order pa- 
rameter TTip on the electron-phonon coupling a for ujq = 1.1, 
K = 0.25, and to = 1. The dashed line is the mean-field result 
and the data with error bar are the Monte Carlo results. 
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